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Abstract. In finance industry portfolio construction deals with how to divide the investors wealth across an 
asset-classes' menu in order to maximize the investors' gain. Main approaches in use at the present are based 
on variations of the classical Markowitz model. However, recent evolutions of the world market showed limi- 
tations of this method and motivated many researchers and practitioners to study alternative methodologies 
to portfolio construction. In this paper we propose one approach to optimal portfolio construction based on 
recent results on stochastic reachability, which overcome some of the limits of current approaches. Given a 
sequence of target sets that the investors would like their portfolio to stay within, the optimal portfolio allo- 
cation is synthesized in order to maximize the joint probability for the portfolio value to fulfill the target sets 
requirements. A case study in the US market is given which shows benefits from the proposed methodology 
in portfolio construction. A comparison with traditional approaches is included. 



1. Introduction 

In finance industry portfolio construction deals with how to divide the investor's wealth across some asset- 
classes selected from a given asset-classes' menu in order to maximize the investor's gain. An asset-class 
is a specific category of investments such as stocks, bonds, cash, and commodities. In quantitative finance 
portfolio construction is achieved by an optimization process. A pioneering work in this regard is the well- 
known Markowitz model [Mar52 j . Given a target performance, Markowitz method provides the optimal 
strategy which minimizes the investment risk. This method and variations of it, are extensively used at the 
present by many asset managers. However, recent evolutions of the world market showed limitations of this 
method and motivated many researchers and practitioners to study alternative methodologies to portfolio 
construction. Main drawbacks of the Markowitz method are summerized in: 

• The asset-classes' performance is assumed to be multivariate gaussian distributed. 

• Investors make a one-shot allocation; the model does not face the portfolio re-balancing during the 
investment life-time. 

The first drawback is shown by basic econometric analysis of the world market to be inappropriate to capture 
relevant statistical properties of asset-classes' performances, while the second one does not allow the investor 
to have a tight control of portfolio evolution during the investment life time. Recently, Optimal Dynamic 
Asset Allocation (ODAA) proposes a methodology to overcome such limitations. ODAA deals with how to 
optimally allocate a multi-period investment. First studies in ODAA faced the problem on how to divide 
the investment among stock and money markets (see e.g. |Mer69j and |Sam69|). More recently, some work 
appeared in the literature concerning optimal strategies for long-lived investors under stochastic investment 
opportunities (see e.g. |BX02j and the references therein). In particular, the work in [BSL97 studies portfolio 
re-balancing in the presence of stochastic variation in the interest rate, the work in [BX02 considers the effects 
of inflation in a portfolio of stock or nominal bonds, the work in [BarOO] and |Xia01] take into account the 
uncertainty in the asset returns prediction. 

In this paper we consider the ODAA problem and we propose a methodology which is based on recent results 
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on stochastic reachability, see e.g. [PLB06 APLS08 Buj04 . Stochastic reachability deals with the synthesis 
of a control strategy which maximizes the probability for the state of a stochastic dynamical control system 
to be in (or to reach) a desired set, representing the set of "good states" . In this paper we reformulate ODAA 
as an appropriate stochastic reachability optimal control problem. Given a specified finite time horizon and 
a sequence of target sets that the investors would like their portfolio to stay within, the optimal portfolio 
allocation is synthesized in order to maximize the joint probability for the portfolio value to fulfill the target 
sets requirements. Main features of this formulation of the ODAA problem are summerized in: 

(Fl) No specific probability distribution assumed on the asset-classes' performances. 

(F2) Sensitivity to market movements. 

(F3) Portfolio evolution control during its life-time. 

We stress that the above features overcome Markowitz limits to portfolio construction. The benefits from 
the proposed approach are illustrated by means of a case study in the US market. A comparison with the 
solution provided by traditional methodologies based on Markowitz optimizer is also included. A preliminary 
investigation on the formalization of the ODAA problem in terms of a stochastic reachability problem appeared 
in the conference publication [PP06J. The present paper provides a detailed and mature description of the 
results announced in RP06 , including econometric analysis of the world market and a case study in the US 
market. 

This paper is organized as follows. In Section [2] traditional approaches to portfolio construction are briefly 
recalled and discussed. In Section [3] an econometric analysis of the main world market is presented, which 
shows limitations of traditional approach to portfolio construction. Section [3] is devoted to the novel approach 
on portfolio construction that we propose. Benefits from this approach are illustrated by means of a case 
study in the US market in Section [5] Details on economic indexes used in such case study are reported in the 
Appendix. Finally Section [6] offers some conclusive remarks and outlook. 



2. An Introduction to Portfolio Construction in Finance Industry 



2.1. Financial Asset— Classes. Portfolio construction deals with how to divide the investor wealth across 
some asset-classes' in order to maximize the investor gain. An asset-class is a specific category of investments 
such as bonds, stocks, and commodities. Assets within the same class generally exhibit similar risk char- 
acteristics, behave similarly in the market-place, and are subject to the same laws and regulations. Due to 
unpredictability of their behaviour, asset -classes dynamics is usually modeled by means of stochastic processes. 
Consider an investment universe made of m asset-classes. Given k £ N define the random vector: 

w k =[w k (l) w k (2) ••• w k (m)] T eR n \ 

whose entries are the performances or equivalently the returns, which are associated with the asset-classes at 
time k. The performance w k (i) is defined as the percentage variation of the asset-class price z k {i) in the time 
interval [k — 1, k], i.e. 

(2- 1 ) w k [i) = — , 

where z k (i) and z k _i(i) correspond to the prices of the i-th asset-class at times k and k — 1, respectively. We 
assume that w k is a correlated random vector. The Expected Return (ER) and Standard Deviation (SD) for 
the i-th asset-class at time k are defined respectively by: 

IM k (i) = E[w k (i)], a k {i) = (E[(w k (i) - ^ k {i)) 2 ]) 1 ' 2 , 

where E[-] is the usual expectation operator. We denote by /i k £ M. m and a k 6 M. m the collection of the 
asset-classes' ERs and SDs at time k, i.e. 

o k = CT fc (l) <7 fe (2) ••• a k (m) f eR m . 



A STOCHASTIC REACHABILITY APPROACH TO PORTFOLIO CONSTRUCTION IN FINANCE INDUSTRY 



3 



Asset-classes' standard deviation is known in finance-industry as volatility and it is a risk measure expressing 
the variability of the asset-class performance around the ER. Covariance matrix (CM) at time k is defined by: 

E fc (i,j) - E[(w k (i) - E[w k {i)]){w k {j) -E[w k (j)})}. 

Quantities ER, SD and CM relate to the first and second moments associated with the random vector w k . 
Univariate higher order statistical indicators of interest are Skewness (SK) and Kurtosis (KU) which are 
defined respectively by: 



e k (i)=E 



w k (i) - fJ. k (i) 
o- k (i) 



and related to the third and fourth moments, respectively. Skewness provides a measure of asymmetry in the 
distribution. Kurtosis quantifies the occurrence of rare events far away from the ER. Gaussian-distributed 
random variable are characterized by a skewness value equal to and by a kurtosis value equal to 3; higher 
values of KU indicate asset-classes with an higher probability of extreme-events, while lower values indicate 
that returns are more predictable and clustered around the ER. Significant deviations from gaussianity are 
usually assessed throughout statistical tests. Given a confidence level CL, which quantifies the accuracy in the 
statistical test, and a data-sample size N = 250, corresponding to the working days per year, an asset-class 
is said to be: 

• Positive Skewed, if the sample skewness SK > X sk ; 

• Gaussian-like Skewed, if the sample skewness SK G]X sk , X sk [; 

• Negative Skewed, if the sample skewness SK < —X sk , 



where: 



X s k — 



6CL 

-l + N 



Moreover an asset-class is said to be: 

• Leptokurtic, if the sample kurtosis KU > 3 + X ku ; 

• Mesokurtic, if the sample kurtosis KU G]3 — X ku , 3 + X ku [; 

• Platykurtic, if the sample kurtosis KU < 3 — X ku , 



where: 



A. 



24 CL 
-l + N 



Each asset-class skewness-kurtosis pair can be represented in the skewness-kurtosis plane. From the above 
taxonomy the skewness-kurtosis plane can be partitioned into nine regions of interest, as depicted in Figure 

HI 



Region 1 
Region 2 
Region 3 
Region 4 
Region 5 
Region 6 
Region 7 
Region 8 
Region 9 



negative skewed and leptokurtic; 
gaussian-like skewed and leptokurtic; 
positive skewed and leptokurtic; 
negative skewed and mesokurtic; 
gaussian-like skewed and mesokurtic; 
positive skewed and mesokurtic; 
negative skewed and platykurtic; 
gaussian-like skewed and platykurtic; 
positive skewed and platykurtic. 
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Figure 1. Skewness-Kurtosis Plane. 



Each of the regions present specific statistical features. For example, Region 1 is characterized by an excess 
of negativ^] rare events, Region 3 by an excess of positive 1 rare events, while Region 5 is not characterized 
by neither a remarkable asymmetry nor by a remarkable fat or thin kurtosis; asset-classes whose skewness 
kurtosis pair are in Region 5 can be assumed to be gaussian distributed. 

2.2. Traditional Approaches to Portfolio Construction in Finance Industry. An asset allocation is 
a vector u E M. m where the i-th entry expresses the amount of the investment in the i-th asset-class. For 
example u — [0.30 0.50 0.20] indicates that we are investing 30% of our wealth in the first asset-class, 50% in 
the second, and 20% in the latter. Let Uk be the asset allocation at time k E N. Some constraints are usually 
imposed on Uk in the investment process. The most common ones in finance industry are: 

• Budget constraint: Y^iLi u k(i) — 1- This equality means that the investor's wealth is fully invested in 
the portfolio; 

• Long-only constraint: Uk(i) > 0, for any i = 1, . . . , m. This inequality implies that short selling is not 
allowed; 

• Risk-budget constraint: (u^-^kUk-i) 1 ^ 2 < o max . This inequality indicates that portfolio risk is 
up-ward bounded by 

We denote by Uk the collection of all asset-allocations allowed in the investment process at time k E N. 
Portfolio performance or equivalently portfolio return r k , can be defined as done for a single asset-class in 



equation (2.1 1, resulting in: 

(2.2) rfc= z fc -z fc -i ; 

Xk-l 

where Xk and Xk-i are the portfolio value at times k and k—1, respectively. It can be shown Luc98 that the 
portfolio performance r^ at time k is given by: 

(2.3) r k = ul^Wf.. 

^Negative and positive events are evaluated with respect to the ER, i.e. an event is negative/positive when its evaluation is 
less/bigger than the ER. 
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The portfolio ER at time k is therefore given by: 

flk = E[rfc] = Ufc_iMfe- 
Analogously the portfolio SD at time k can be expressed by: 

& k = (E[(r fe - p k ) 2 ]fl 2 = (uL^fc-i) 172 . 

Traditional approaches to portfolio construction in finance industry make use of (variations of) the classical 
Markowitz method [M ar52] . Markowitz method provides an optimal solution to the asset-classes' allocation 
problem. The optimality criterion considered is the risk minimization. Let us consider an asset -classes' menu, 
characterized by ERs CM and a target return r. Investors in Markowitz world wish to minimize 
portfolio risk among the portfolios ensemble with target return r. The Markowitz optimization problem is 
therefore formalized as follows: 

min^gc/fc ul T, k u kl 



(2.4) 

Feasible target returns f belong to the interval [f m i n ,r max ], where: 



= (arg inf u k H k u k ) p, k , r max = sup u eU .u%p k . 

u. k eu k 



By performing the optimization problem in ( |2.4[ ), with target return r ranging in [r OT j n , r max ], we obtain 
the collection of all optimal portfolios, known in the literature as the Efficient Frontier [Lue98 . Efficient 



algorithms are known in the literature for solving the optimization problem in (2.4 1. While being handablc 



for practical computation, Markowitz method exhibits some drawbacks, summarized in: 

• The asset-classes' performance is assumed to be multivariate gaussian distributed: the only risk-figures 
involved in the portfolio optimization are contained in the covariance matrix CM. 

• Investors make a one-shot allocation, the model does not face the problem of how (optimally) rebalance 
portfolio during the investment life- time. 

In the next section we will show the impact of such drawbacks in portfolio construction. 



3. Limits of Gaussian Models in Portfolio Construction 

We start by providing an academic example which illustrates implications of gaussian-based asset-classes 
modeling in portfolio construction. Consider an asset-classes' ensemble composed of three un-correlated 
synthetic asset-classes. The asset-classes returns are modeled as follows: 

(3.1) Wk(l) = 7) u>/c(2) = 2p - 7, w k (3) = r), 

where p € R, random variable 7 is r(a, 0) distributed with a = 1 and f3 = p, and 77 is a gaussian random 
variable with ER and SD equal to p. In this example we set p = 0.03. Figure [2] illustrates the marginal 



probability density functions of the random variables in (3.1 1. 
It is possible to show that the considered asset-classes share the same ER and SD, although they are very 
different in terms of skewness and kurtosis. The following table reports statistical indicators associated with 
w k (l), w k (2) and w k (3): 



(3.2) 



Asset 


ER 


SD 


SK 


KU 


Asset 1 


0.03 


0.03 


2 


9 


Asset 2 


0.03 


0.03 


-2 


9 


Asset 3 


0.03 


0.03 





3 



With reference to the skewness-kurtosis plane in Figure [T] Asset 1 is in Region 3, Asset 2 is in Region 1 
and Asset 3 in Region 5. Since ER and SD of the asset-classes are the same, asset allocation based on these 
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-0.2 -0.15 -0.1 -0.05 



FIGURE 2. Synthetic data-set. Marginal probability density functions. 



statistical indicators is not able to distinguish among asset-classes. Indeed, Markowitz optimizer selects the 



optimal solution Uq of (2.4 1 given by the equally-weighted portfolio, i.e. 



1/3 1/3 1/3 



On the other hand a simple inspection of the table in (3.2 1 indicates that the three asset-classes present 



important differences in the marginal probability distribution. A prudent investor would invest in Asset 1 
which is characterized by low probability to suffer negative returns but low probability to reach large positive 
returns. On the way around Asset 2 is much riskier than Asset 1 but it guarantees higher probability to 
get large positive returns. This simple example shows limitations of portfolio construction based on gaussian 
modeling of asset-classes. On the other hand deviations from gaussianity are evident in the main world market 
asset-classes. 

In the following we report a basic econometric analysis of the main world markets. We analyze 33 indices 
from January 1-st 1992 to December 31-th 2008. Indices are in local currencies (see Appendix A) and daily- 
based (closing price values) . In order to test the gaussian hypothesis we perform the Jarque-Bera test to 95% 
Confidence Level on the indices' returns. We considered 1-year non-overlapping windows from 1992 to 2008; 
this choice allows us to monitor the time-evolution of higher order statistics year by year. 
Results are shown in Figure [4] Asset-classes in Region R5, for which gaussinity hypothesis can be assumed 
account at a maximum of 33% in 1999 and a minimum of 0% in 2008. Asset-classes in Regions Rl and R3 
account at a maximum of 72% in 1996 and a minimum of 27% in 1999, 2002 and 2005. Some meaningful 
specific cases are reported below: 



JPY-EUR. Figure 3(a) depicts statistical historical data associated with JPY-EUR in the skewness- 
kurtosis plane. Each dot represent a skewness-kurtosis pair in each year ranging from 1992 to 2008. 
Only 5 over 17 dots are located in Region 5; 

Emerging Market Govy. Figure [3 (b) | depicts statistical historical data associated with Emerging Mar- 
ket Govy in the skewness-kurtosis plane. Each dot represent a skewness-kurtosis pair in each year 
ranging from 1992 to 2008. Only 1 over 17 dots is located in Region 5. 



This econometric analysis highlights the evidence and the impact of higher order statistics in common-used 
risk measure. 
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(a) JPY-EUR in the skewness-kurtosis plane. 



(b) Emerging Market Govy in the skewness-kurtosis plane. 



Figure 3. JPY-EUR and Emerging Market Govy in the skewness-kurtosis plane. Each dot 
represents a skewness-kurtosis pair in each year ranging from 1992 to 2008. 

4. A Novel Approach to Portfolio Construction 



In this section we illustrate an alternative approach to portfolio construction, based on stochastic optimal 
control. A stochastic dynamical control system modeling the portfolio value dynamics can be derived by 



combining equations in (2.2 1 and (2.3 1, resulting in: 

(4.1) x k +i = x k (l + v%w k+ i), 

where: 



k e N, 



• Xk G X = K is the state, representing the portfolio value at time k; 

• Uk E Uk C E m is the control input, representing the portfolio allocation at time k; 

• Wk £ K m is a random vector describing the asset-classes' returns at time k E N. 



Let (fl,F,P) be the probability space associated with the stochastic system in (4.1|. Portfolio value Xk at 
time k = is assumed to be known and set to xq — 1. The mathematical model in (4.1 1 is characterized by: 



(Fl) No specific distribution on the asset-classes' returns. 

As already discussed in the previous section this feature is important in asset-allocation. We model asset- 
classes' returns by means of Mixture of Multivariate Gaussian Models (MMGMs), which provide an accurate 
modeling of non-gaussian distributions while being handablc for practical implementations. We recall that a 
random vector Y is said to be distributed according to a MMGM if its probability density function py can 
be expressed as the convex combination of probability density functions py i of some multivariate gaussian 
random variables Y;, i.e. 



N 

py{v) = ^2\PYi(y), 
i=i 



N 



A* €[0,1], 2^ = 1. 



Some further constraints are usually imposed on coefficients Xi so that the resulting random variable Y is well 
behaved, by requiring for example semi-definiteness of covariance matrix and unimodality in the distribution. 
The interested reader can refer to |BSS08j for a comprehensive exposition of main properties of MMGMs. 
The class of control inputs that we consider in this paper is the one of Markov policies |Ber01j . Given a finite 
time horizon iVgNa Markov policy is defined by the sequence 

7T = {U ,U 2 , ...,U N -l} 
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Figure 4. Econometric analysis of the main world market. 



of measurable maps Uk '■ X Uk- Denote by Uk the set of measurable maps Uk '■ X — > Uk and by TIn 
the collection of Markov policies. For further purposes let be ir k — {uk, Wfc+i, itjv-i}- Markov policies 
adequately model portfolio allocations, since they allow the investor to have: 

(F2) Sensitivity to market movements. 

Indeed, control input Uk depends on Xk and hence on the market performance at time t € [k — 1, k]. 

Let us consider a finite time horizon TV which represents the life-time of the considered investment. Our 
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approach in the portfolio construction deals with how to select a Markov policy 7r in order to fulfill some 
specifications on the portfolio value x k at times k = 1, . . . , N. The specifications are defined by means of a 
sequence of target sets {Xi, X%, X^} with Xi C X. The investor wishes to have a portfolio value Xk at time 
k that is in X%. Typical target sets Xk are of the form Xk — [x fe , +00 [ and ask for achieving a performance 



that is downside bounded by x fc € K. An example of target sets pattern is illustrated in equation (5.1 1 and 
models an investor specification which requires to beat a target return 7% (annualized value) at maturity at 
the end of the investment life time. This formulation of specifications allows the investor to have a 

(F3) Portfolio evolution control during its life-time, 

since target sets Xk depend on time k. The portfolio construction problem is then formalized as follows: 

Problem 4.1. (Optimal Dynamic Asset Allocation (ODAA)) Given a finite time horizon N G N and a 
sequence of target sets 

(4.2) {X 1 ,X 2 ,...,X N }, 

where Xk are Borel subsets of X , find the optimal Markov policy it that maximizes the joint probability quantity 

(4.3) p({uj e ft ■. x e x ,xi e x 1: ...,x N e x N }). 



The ODAA Problem can be solved by using a dynamic programming approach [BcrOl and in particular 
by resorting to recent results on stochastic reachability (see e.g. |PLB06[ |APLS08, !Aba071 IBL031 |Buj04| ). 
Since the solution of Problem |4.1| can be obtained by a direct application of the results in the work of 
PLB06 APLS08 in the following we only report the basic facts which lead to the synthesis of the optimal 
portfolio allocation. Given x € X and u € K m denote by Pf( x .u,w k ) the probability density function of random 
variable: 

f(x, u, w k+1 ) = x(l + u T w k+1 ), 
associated with the dynamics of the system in ( |4.1[ ). Given the sequence of target sets in ( |4.2| and a Markov 
policy 7r we introduce the following cost function V which associates a real number V(k, x, ir k ) e [0, 1] to a 
triple (k, x, ir k ) by: 

vf . fc , _ f Ix h (x), if k = N, 

V{>c,x,Tr f x ^ V ( k + 1>Z)7r k+i) Pnxumi) ( z ) dz , if k = N - 1, N - 2, 1, 0, 

where Ix N (x) is the indicator function of the Borel set Xn, i.e. Ix N (x) = 1 if x G Xn and Ix N {x) = 0, 
otherwise. Results in [PLB06 show that cost function V is related to the probability quantity in (4.3 1 as 
follows: 

p({lj eSl:x e4iie X 1: ...,x N e x N }) = V{Q,x ,tt). 

Hence the ODAA Problem can be reformulated, as follows: 

Problem 4.2. (Optimal Dynamic Asset Allocation) Given a finite time horizon N S N and the sequence of 
target sets in (4-2) compute: 

7r* = arg sup V(0,Xo,Tt). 

Tt£Yl N 

The above formulation of the ODAA Problem is an intermediate step towards the solution of the optimal 
control problem under study which can now be reported hereafter. 

Theorem 4.3. [PLB06] The optimal value of the ODAA Problem is equal to 

P* = Jq( x o), 

where Jq{x) is given by the last step of the following algorithm, 
j n{x) = I Xn (x), 

( ' ' Jk{x) = sup UkEUk J Xk+i J k+1 (z)p f{XiUkiWk+l) {z)dz, k = N-l,N-2, ...,1,0. 
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C 


B 


E 


ER (ami) 


3.24% 


5.46% 


10.62% 


SD (ann) 


0% 


4.45% 


14.77% 


SK 





-0.46 


-0.34 


KU 


3 


4.25 


5.51 



correlation 


C 


B 


E 


C 


1 








B 





1 


0.0342 


E 





0.0342 


1 



Table 1. Statistical indicators of the asset-classes menu. 



The above result proposes an algorithm whose outcome is precisely the optimal Markov policy solving the 
ODAA Problem. A financial strategy interpretation of the above result can be found in the so-called contrarian 
strategy, see e.g. |Cha88j and the references therein. A financial strategy in asset allocation is said to be 
contrarian when it suggests to buy risk asset-classes when the market is down and to sell them when the 
market is performing well. In the next section we will illustrate the contrarian attitude of the proposed 
approach in a case study. 

We conclude this section by briefly discussing computational issues related with the algorithm proposed in 



Theorem 4.3 Although as discussed in Aba07], algorithms for studying stochastic reachability for general 
stochastic nonlinear (and hybrid) control systems are demanding from the computational point of view, the 
computational complexity of the proposed algorithm remains tractable because of the following reasons: (i) the 
state variable of the system in ( |4.1| is scalar; (ii) properties of MMGM, that we use in modeling distribution of 
Wk, allow computational efficient solutions in solving the basic optimization problem in ( |4.4[ ) (see e.g. |BSS08j ). 
Moreover we stress that the optimization process required in Theorem |4.3| to portfolio construction can be 
done off-line and that design of financial products does not require tipically, severe time constraints. 



5. A Case study: A Total Return Portfolio in the US Market 

In this section we illustrate the proposed methodology by designing an asset allocation in the US market. 
Let us consider a total-return US fund-manager. The investment's universe consists of three asset-classes: 
Money-market, US Bond, and US Equity markets. Details on the indices used in the analysis arc reported in 
the table below: 



Label 


Asset 


Index 


C 


Money market 


US Generic T-bills 3 months 


B 


US Bond 


JP Morgan US Government Bond All Maturity 


E 


US Equity 


S&P500 



Money market, bond and equity are tipically used as a portfolio ensemble in asset -allocations because of their 
diversified statistical properties which allow the investors to deal with different financial scenarios. 
Time-series are in local currency (US dollars) and weekly-based from January 1-st 1988 to December 28-th 
2007. Statistical indicator^Jin the selected portfolio menu are reported in Table[l] By comparing asset-classes 
ERs it is readily seen that US Equity (asset-class E) ensures higher performances than US Bond (asset-class 
B) which in turn, ensures higher performances than Money Market (asset-class C). On the other hand, asset- 
classes SDs show that US Equity is riskier than US Bond which is riskier than Money Market. Moreover, 
skewness and kurtosis reported in TableJTJshow that asset-classes present significant deviations-to-gaussianity 
(Jarque-Bera test; 99% Confidence Level). Bond and Equity markets are leptokurtic and negative skewed; 
with reference to Figure [l] while Asset C is in Region 5, Assets B and E are located in Region 1. 
The asset allocation that we want to synthesize is characterized by a 2-years life-time. We consider a weekly 
re-balancing, i.e. we suppose to re-balance our portfolio once per week. The finite time horizon associated 
with this asset allocation is therefore given by N = 104. The investor specification is to beat a target return 



Standard deviation, skewness, kurtosis and correlation have been estimated by considering the historical average in the 
selected time-window daily observation. Expected return has been computed assuming a constant Sharpe-ratio ILue 98 of 0.50. 
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7% (annualized value) at maturity; his risk-budget corresponds to 7% (ex-ante) monthly Value-at-RislJ^] at 
99% Confidence-Level. 

The above specification translates in constraints on the set Uk of portfolio allocations and on an appropriate 
choice of the sequence of target sets X%. With regard to the definition of Uk, a portfolio allocation Uk £ Uk is 
required to satisfy budget, long-only and risk-budget constraints with a max = VaR Vl~2/2.3263 = 0.1042, see 
footnote ( 3 ). Moreover the target sets Xs' formalization results in: 

X Q = {!}, 

X k = [0,+oo[, fc=l,...,103, 
(5.1) X Wi = [1.07 2 ,+oo[. 

The optimization criterion consists in maximizing the probability 

P({uj e : x 10i > 1.07 2 }), 
where x W4: is the portfolio value at the end of the second year (investment life-time). 

The first step in the methodology illustrated in the above section consists in finding a MMGM which ap- 
propriately describes the market scenario behaviour. The chosen MMGM Y is composed of two multivariate 
gaussian random variables Y\ and Y 2 , characterized by the following univariate statistics: 

iT 



0.000611 
0.000069 



0.001373 
0.005666 



0.002340 
0.019121 



(5.2) 

and same correlation matrix: 



0.000683 -0.016109 -0.017507 
0.000062 0.006168 0.052513 ] T . 





C 


B 


E 


c 


1 


0.0633 


0.0207 


B 


0.0633 


1 


-0.0236 


E 


0.0207 


-0.0236 


1 



Symbols /j, 1 and cr l indicate (resp.) the ER and SD of the z-th multivariate gaussian model, which are associated 
with the three asset-classes. Values are weekly-based. The MMGM parameters expressing the probability 
related to each gaussian world are: 



Ai = 0. 



A 2 = 0.02. 



The obtained MMGM correctly represents the univariate statistics of the asset-classes up to fourth-order and 
the correlation patterns up to the second-order. The above MMGM has been obtained by an optimization 
process which minimizes the errors between the univariate four moments associated with the time series of the 
asset-classes involved and the selected MMGM; the obtained error resulted in 2.24E 1 — 7. Additional constraints 
of semi-dcfinitcncss in the covariance matrix and unimodality of the univariate marginal distributions have 
been considered in the selection of the MMGM. 

By applying the methodology illustrated in the previous section the optimal portfolio allocation is synthesized. 
The optimal Markov policy at time k = 0, corresponding to the initial condition xq = 1, is given by 29.50% 
Money market, 0% Bond market, and 70.50% Equity market, as shown in the first panel of Figure [5] After 
the first week fund-manager revises the portfolio allocation. Optimal Markov policies corresponding to the 
end of each semester and 1-week before the investment maturity (k = 26, 52, 78, 103) are reported in Figure|5] 
In the second panel abscissas report the portfolio value x 2 6 at time k = 26. For each portfolio realization x 2 6i 
the map gives the corresponding portfolio allocation. For example a portfolio value of x 2 6 — 1-000 indicates 



^Value-at-Risk (VaR) quantifies the potential loss of a portfolio to a given time horizon and to a certain confidence level. 
Parametric evaluation of the VaR permits to relate it to the portfolio volatility | Lue98| ; for example VaR to 1 month horizon and 
99% confidence level is 2.3263 SD /Vl2, where SD indicates the portfolio standard-deviation. 
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Figure 5. Optimal portfolio allocation at times k = 1, 2, 8. 



to allocate about 30% of the investors' wealth in US Bond (asset-class B) and 70% in US Equity (asset-class 
E). It is readily seen that as the portfolio value decreases the portfolio allocation gets riskier and riskier; 
for example optimal policy yields 0% Cash for portfolio value less than 1.0081. This means that until this 
threshold is reached, the model needs to allocate the maximum allowed budget risk expressed by the Value 
at Risk 99% constraint. In the other way around as performances gets better and better the model suggests 
to reduce the risk exposure: for a portfolio value greater than 1.0919, a 100% Cash allocation guarantees to 
reach the performance target. The attitude of this optimal strategy to increase risk exposure in presence of 
portfolio draw-downs, and conversely to reduce it in case of positive performance is known in the literature 



A STOCHASTIC REACHABILITY APPROACH TO PORTFOLIO CONSTRUCTION IN FINANCE INDUSTRY 



13 



0.1 
0.09 
0.08 
0.07 
0.06 



O 0.05 
o. 



0.04 
0.03 
0.02 
0.01 



0^ 

-0.1 



-0.05 



0.05 



0.1 



0.15 



0.2 



0.25 



0.3 



FIGURE 6. Return Localization. 

as contrarian strategy [Cha88 . Optimal Markov policies at times k = 52,78 exhibit similar characteristics to 
the one for k = 26, as shown in Figure [5| It is readily seen that as k increases the portfolio re-balancing gets 
sharper and sharper: as time passes by, portfolio allocation strategy becomes more and more aggressive and 
indicates to buy more and more US Equity (asset-classes E) in the attempt to reach the target goal. Main 
differences is that the optimal Markov policy require to be riskier than the one corresponding to the previous 
map. Let us consider a portfolio value equal to 1.06. At k — 26 the optimal allocation is 62% Cash, 12% 
Govy, 26% Equity. After a semester (k = 52), the optimal allocation for the same portfolio value is 16% Cash, 
26% Govy, 58% Equity. At k = 78 the optimal allocation is 0% Cash, 30% Govy, 70% Equity. The last map 
at k = 103 presents similar feature. It should be stressed that the optimal solution is particularly sharp. For 
a portoflio value equal to 1, 1435 the model suggests to use all the available budget risk, while for a portfolio 
value equal to 1, 1446 a full Cash allocation ensures to reach the target. 

The maximal probability of achieving the investment goal is: 

p* = J a {x ) = 77.62%. 

In order to make a validation of this result we run a Montecarlo simulation with 10 6 scenarios. The results 
yield a probability of 77.76% and are depicted in Figure [6] The different probability values achieved by the 
optimization process and the Montecarlo simulations are due to the numerical errors in the optimization process 
and to the number of realizations chosen in the Montecarlo simulations; nevertheless the two probability values 
are close one each other. Figure [6] shows that most of the reali zatio ns are placed at the right of the portfolio 



value x = 1.07 2 as required by the target sets specifications in (4.2 1. 



In order to make a comparison with traditional standard approaches in asset-allocation, we run the same 
exercise for a Markowitz investor with a finite time horizon TV = 104. In this case the optimal solution 
requires to invest the full budget-risk, corresponding to an allocation of 0% Cash, 30% Govy, 70% Equity. 
This allocation has been derived at first by optimizing according to the mean-variance paradigm, and then 
by selecting the optimal portfolio on the efficient frontier which maximizes the probability to beat the target 
return in two years horizon. The resulting probability of achieving the investment target goal is 61.90%. This 
probability has been evaluated by performing a Montecarlo simulation with 10 6 scenarios. By comparing the 
probability of achieving the target goal in the two approaches it is readily seen that the ODAA approach 
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provides a differential of f5.86% with respect to the Markowitz approach, which is remarkable in portfolio 
allocation in finance industry. 



6. Conclusion 



In this paper we addressed Optimal Dynamic Asset Allocation. Given a specified finite time horizon and 
a sequence of target sets that the investors would like their portfolio to stay within, the optimal portfolio 
allocation is synthesized in order to maximize the joint probability for the portfolio value to fulfill the target 
sets requirements. The proposed approach has been shown to overcome asset-classes gaussian-based modeling 
and static portfolio allocation limits of the methodologies currently in use in finance industry. 
Contrarian strategy attitude of the proposed methodology can result in being rather aggressive if applied to 
some specific financial products in case of large market draw downs. The authors are currently investigating 
some approaches based on stochastic hybrid systems modeling |Buj04| with the goal of making the proposed 
methodology suitable of application to a wider range of financial products. 
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Appendix: time series-providers 



Data-providers and details of the time-series listed in Figure [4] are: 

• FX rates. Bloomberg; 

• Government Bonds. JP Morgan (All Maturity indices); 

• Inflation-Linked Bonds. Barclays Capital International (All Maturity and Total Return indices); 

• Corporate Bonds and High Yields. Merril Lynch; 

• Equity markets. Morgan Stanley Capital International; 

• Convertible. UBS (Total Return and At-The-Money indices); 

• Real Estate. EPRA/NAREIT (Total Return and At-The-Money indices); 

• Commodities. Goldman Sachs (Total Return and At-The-Money indices). 

All time-series are in local currencies, Commodities are in USD. Bloomberg |BLO] code of the asset-classes 
listed in Figure [4] are (resp.): 



gbpeur curncy, 
jpmgemlc index, 
bciult index, 
msdlemu index, 
msdlsg index, 
mxemsc index, 
mclasg index, 
ugna index, 



jpyeur curncy, 
jpmtus index, 
cOaO index, 
msdlus index, 
msclcgfm index, 
mcldus index, 
ucbiemae index, 
efjp index, 



usdeur curncy, 
jnycjp index, 
crOO index, 
msdluk index, 
mselegfl index, 
mclduk index, 
ucbijvaj index, 
gscitrsi index. 



cadeur curncy, 
bceelt index, 
jOaO index, 
msdljn index, 
mselegfa index, 
mclajn index, 
ucbiusam index, 



audeur curncy, 
bcitlt index, 
heOO index, 
msdlhk index, 
mseltcf index, 
mclahk index, 
exuk index, 
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